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dimensional Abelian ( Non-Abelian ) Chern-Simon theory 



O 

o 

(N 



Jinwu Ye 

Department of Physics and Material Research Institute, The Pennsylvania State University, University Park, PA, 16802 

(February 1, 2008) 

By applying the simple and effective method developed to study the the gauge-invariant fermion 
Green function in 2 + 1 dimensional non-compact QED, we study the gauge- invariant Green function 
in 3 + 1 dimensional QED and 2 + 1 dimensional non-compact Chern-Simon theory. We also extend 
our results to the corresponding SU(M) non-Abelian gauge theories. Implications for Fractional 
Quantum Hall effect are briefly discussed. 
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I. INTRODUCTION 

In any systems of gauge fields ( Abelian or Non-Abelian 
) coupled to matters ( fermions or bosons ), the conven- 
tional Green function is defined as: 



G(x,y) =< ip(x)il>(y) > 



(1) 



In momentum space, the fully-interacting Fermion 
Green function in Eqn.l takes the form G(k) = T^pf ■ 
In 4 ( 3 ) dimensional space-time, it corresponds to JiiHa 
( ) where the anomalous dimension rj can be calcu- 
lated by standard Renormalization Group (RG) by ex- 
tracting UV divergences. Unfortunately, this Green func- 
tion is not gauge invariant, r\ depends on the fixed gauge 
in which the calculation is done. 

Schwinger proposed the following Gauge-invariant 
Green function [1]: 



G' 



'(x,y) =< ip(x)e ie J* af > i6)d ^Tp(y) > 



(2) 



where the inserted Dirac string makes the Schwinger 
Green function gauge invariant. Q mv depends on the 
the integral path C from x to y. For simplicity reason, 
we take C to be a straight line [2] . 

In fact, the path C should be determined by the un- 
derlying physical systems. In QED or QCD, in principle, 
one should perform average over all possible paths from 
x to y with some weights, but what kind of weight factors 
should be used is still unclear. Non-smooth paths with 
cusps or intersections will also cause additional complica- 
tions. Perturbative QCD makes sense only at short dis- 
tance, so choosing straight line may be reasonable. See 
[3] for some preliminary discussions. In condensed mat- 
ter systems such as Fractional Quantum Hall Effects or 
High Temperature Superconductors (HTS), formulated 
on a lattice, the physical quantities such as real elec- 
tron Green function in FQHE or Angle Resolved Photo- 
Emissions Spectroscopy (ARPES) in HTS are indepen- 
dent of path, therefore taking a straight line is not only 
simplest and also plausible. In [29], we will discuss the 
path dependence of gauge-invariant Green functions in 
different condensed matter systems formulated on lat- 
tices. 



In Quantum Chromo dynamics (QCD), the Schwinger 
gauge-invariant Green function is closely associated with 
the hadronic bound states and was studied before [3] . In 
recent years, its importance in condensed matter system 
was recognized in the contexts of Fractional Quantum 
Hall systems [4] and high temperature superconductors 
[5-7]. 

By a singular gauge transformation which attaches two 
flux quanta to each electron, an electron in an exter- 
nal magnetic field was mapped to a composite fermion 
in a reduced magnetic field [8-10]. Although trans- 
port properties which are directly related to two particle 
Green functions can be directly studied in the composite 
fermion language, the tunneling density of states which is 
directly related to the single particle electron Green func- 
tion is much more difficult to study. In fact, the single 
particle electron Green function is equal to the gauge- 
invariant Green function of the composite fermion which 
was evaluated for non-relativistic fermions by phenomc- 
logical arguments in Ref. [4]. 

Most recently, the importance of gauge-invariant 
Green function of fermion to Angle Resolved Photo- 
Emission (ARPES) data in high temperature supercon- 
ductors was discovered independently by Rantner and 
Wen (RW) [5] and the author [6] in different contexts. 
Starting from U (1) or SU (2) gauge theory of doped Mott 
insulators [11-13], RW discussed the relevance of this 
gauge-invariant Green function to ARPES data. They 
also pointed out that in temporal gauge, the equal-space 
gauge-invariant Green function in Eqn.2 is equal to that 
of the conventional gauge dependent one in Eqn.l. Start- 
ing from a complementary (or dual ) approach pioneered 
by Balents et al [14], the author studied how quantum [7] 
or thermal [6] fluctuations generated hc/2e vortices can 
destroy c£-wave superconductivity and evolve the system 
into underdoped regime at T — or pseudo-gap regime 
at finite T. In the vortex plasma regime around the fi- 
nite temperature Kosterlitz-Thouless transition [15], the 
vortices can be treated by classical hydrodynamics. By 
Anderson singular gauge transformation which attaches 
the flux from the classical vortex to the quasi-particles 
of d-wave superconductors [16-19], the quasi-particles 
(spinons) are found to move in a random static magnetic 
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field generated by the classical vortex plasma [6]. The 
electron spectral function G(x,t) =< C a (0,0)C^(x,t) > 
is the product of the classical vortex correlation function 
and the gauge invariant Green function of the spinon 
in the random magnetic field. Tcchniqually, this static 
gauge invariant Green function is different from the orig- 
inal dynamic Schwingcr gauge invariant Green function. 
Conceptually, both arc single particle gauge-invariant 
Green function and are physically measurable quantities. 

Recently, two different methods are developed to cal- 
culate the gauge invariant Green function Eqn.2. The au- 
thor in Ref. [20] evaluated it in a path integral representa- 
tion. In Ref. [21], by applying the methods developed to 
study clean [22] and disordered [23] FQH transitions and 
superconductors to insulator transitions [24], the author 
developed a very simple and effective method to study the 
gauge invariant Green function. In the context of 2+1 di- 
mensional massless Quantum Electro-Dynamics (QED3) 
[25,12,13], the author studied the gauge-dependent Green 
function both in temporal gauge and in Coulomb gauge. 
In temporal gauge, the infrared divergence was found 
to be in the middle of the contour integral along the 
real axis, it was regularized by deforming the contour 
into complex plane by physical prescription, the anoma- 
lous dimension was found and argued to be the same 
as that of gauge-invariant Green function. However, in 
Coulomb gauge, the infra-red divergence was found to 
be at the two ends of the contour integral along the real 
axis, therefore un-rcgulariablc, anomalous dimension was 
even not defined. This infra-red divergence was shown 
to be canceled in any physical gauge invariant quanti- 
ties such as p function and correlation length exponent 
v as observed in Refs. [22-24]. The author also studied 
the gauge-invariant Green function directly by Lorentz 
covariant calculation with different gauge-fixing param- 
eters and find that the exponent is independent of the 
gauge fixing parameters and is exactly the same as that 
found in temporal gauge. 

In this paper, we apply the simple and effective method 
developed in Ref. [21] to study the two interesting phys- 
ical systems: 3 + 1 dimensional QED and 2+1 di- 
mensional Chern-Simon theory and their corresponding 
SU(M) non-Abelian counter-parts. The importance of 
the first system is obvious in high energy physics. The 
second system is closely related to the high temperature 
behaviors of 3 + 1 dimensional QED or QCD. It may 
also describe Fractional Quantum Hall (FQH) transitions 
[10,27,22,23]. The gauge invariant Green function Eqn.2 
is a relativistic analog of tunneling density of state in 
FQH system studied by phenomelogical method in [4]. 

The paper is organized as follows. In the next section, 
we study QED4 in both temporal gauge and Lorentz co- 
variant gauge, we also extend our results to non-Abelian 
SU(M) QCD. In section III, we study QED3 with Chern- 
Simon term also in both temporal gauge and Lorentz co- 
variant gauge, we also extend our results to non-Abelian 



SU(M) Chern-Simon theory. Finally, we reach conclu- 
sions by summarizing our results in three simple and in- 
tuitive rules. 

II. 3 + 1 DIMENSIONAL QED 

The standard 3+1 dimensional massless Quantum 
Electro- Dynamics ( QED4 ) Lagrangian is: 

£ = -ipalnidn ~ iea^a + ^(f^) 2 (3) 

Where tp a is a four component spinor, a = 1, • • • , N are N 
species of Dirac fermion, 7^ are 4x4 matrices satisfying 
the Clifford algebra {7^,7^} = —28^. 

In contrast to QED3 studied in [21], the coupling con- 
stant e is marginal ( or dimensionless) at 4 space-time 
dimension. Straight-forward perturbation suffices. 

With the gauge fixing term ^-(c^a^) 2 , the gauge field 
propagator is: 

D» v = ^{5» v -{l-a) k -^) (4) 

By extracting UV divergence, we find the anomalous 
exponent for the gauge dependent Green function in 
Eqn.l 



Note that in Landau gauge a — 0, r\ vanishes ! In 

usual textbooks, the detailed calculations were given in 
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Feymann gauge a = 1 with r/ = — Obviously, r\ 
is a gauge-dependent quantity and its physical meaning 
in any Lorentz covariant gauges is not evident. In the 
following section, we will calculate r\ in two Lorentz non- 
covariant gauges: temporal gauge and Coulomb gauge. 

A. The calculation in Temporal gauge 

As stated in the introduction, the main focus of this 
paper is the gauge-invariant Green function Eqn.2. In 
temporal gauge, the equal-space gauge invariant Green 
function is the same as the gauge dependent one [5]. 

The strategy taken in Ref. [21] is to calculate the con- 
ventional Green function in temporal gauge = and 
then see what we can say about the gauge invariant Green 
function. As shown in [21] in the context of QED3, in 
temporal gauge and Coulomb gauge which break Lorentz 
invariance, we run into both UV and IR divergences. A 
sensible and physical method was developed to regular- 
ized these plaguey IR divergences. Here we take the same 
strategy and apply similar method to regularized these 
IR divergences in the context of QED4. 
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With the notation K = (fco, k), in do = gauge, it is 
easy to find the propagator: 



(6) 



The one-loop fermion self-energy Feymann diagram is 




Fig 1: The fermion self-energy diagram in temporal gauge 

The corresponding expression is 

d 3 Q ^(K-Q)^ 1 



(2tt)3 74 (K-Q) 2 U Q 



1 2 o } 



(7) 



By using standard 7 matrices Clifford algebra and sup- 
pressing the prcfactor — «e 2 , we can simplify the above 
equation to: 



E(tf) = 70 



+ 7, 



d 4 Q (k - q) 



:(3- 



2 / 



(27r)4(^-0) 2 g 2 ^ 9o 
d 4 Q (k + g)j 
(2tt) 4 (K - Q) 2 Q 2 
d 4 Q q 2 (k + q)j - 2q- kq { 
(2tt)4 (K - Q) 2 Q 2 ql 



(8) 



We choose the external momentum K to be along 
z axis, then Q4 = Qcos9, Q3 = Q cos 9\ sin#, Q2 = 
Q cos </> sin 0i sin^, Qi = Q sin sin 6*i sin0, d A Q = 
Q 2 dQd<j) sin 9\d9i sin 2 9d8. Setting x = — cos 9, we find 
the logarithmic divergence: 



7o fc o 
4^3 



^•1 

J dx V 1 - x 2 (-Ax 2 + x- 2 ) log A 



(9) 



The integral can be rewritten as: 



dx^- 



7T 

2 



(10) 



As expected, we run into IR divergence at x — which 
is in the middle point of the contour integral on the real 
axis from —1 to 1. Fortunately, by physical prescription, 
we can avoid the IR singularity at x = by deforming 
the contour as shown in Fig. 2 



Fig 2: The contour path C to bypass the IR singularity in tem- 
poral gauge 



The divergent part in Eqn.10 becomes: 
a/! - 1c 2 



/_", 



dz- 



x* 



VT^ 2 



(ii) 



Putting back the prefactor 
swer: 



we get the final an- 



. 2 7 fc 3?r 3e 2 

~ ie l^ ( ~T )logA = J7ofc °^2 logA 

We can identify the anomalous dimension as 



r\ = 



3e2_ 
8ir 2 



(12) 



(13) 



We expect this is the correct anomalous dimension of 
the gauge-invariant Green function Eqn.2 in QED4. It 
is consistent with the result in [3,20] achieved with very 
different method. 

In the above calculation, we choose a cut-off A in 4- 
momentum Q, Just as in QED3 studied in [21], we can 
introduce an alternate cut-off A only in 3-momcntum q, 
but integrate the frequency go freely from —00 to 00. Us- 
ing the similar infra-red regulization scheme as in Fig. 2, 
we find exactly the same answer as Eqn.13. This agree- 
ment indeed shows that the exponent Eqn.13 is universal 
independent of different cut-offs or different renormaliza- 
tion schemes. 

Just as in QED3 studied in [21], we could evaluate the 
Green function in Coulomb gauge diai = 0. As in tempo- 
ral gauge, we also run into IR divergence. In both cut-off 
Q < A and cut-off q < A, we run into IR divergences at 
x = ±1 which are at the two end points of the contour 
integral on the real axis from —1 to 1. Unfortunately, 
from physical prescription, we are unable to avoid the IR 
singularities at the two end points x — ±1 by deform- 
ing the contour. Therefore, we arc unable to identify 
anomalous exponent. Furthermore, the results are dif- 
ferent in the two different cutoffs. This should cause no 
disturbance, because, in contrast to the Green function 
in temporal gauge, the Green function Eqn. 1 in Coulomb 
gauge does not correspond to any physical quantities. 
All these IR divergences should disappear in any gauge- 
invariant physical quantities like j3 function and critical 
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exponents. The final answers for these physical quan- 
tities should also be independent of different cut-offs or 
different rcnormalization schemes. 



B. Lorentz covariant calculation 

In this subsection, we will calculate the gauge invari- 
ant Green function directly in Lorentz covariant gauge 
Eqn.4 without resorting to the gauge dependent Green 
function. We will also compare with the result achieved 
in temporal gauge. 

The inserted Dirac string in Eqn. 2 can be written as: 



(14) 



where the source current is: 



j» = j c dr^6(x,-Ur)) (15) 

where r parameterizes the integral path C from x to y. 

We will follow the procedures outlined in detail in 
Ref. [21]. (1) Combining the source current with the 
fermion current j^(x) — tp(x)'-f fl ^p(x) to form the total 
current f^x) — j^(x) + j s Ax) (2) Integrating out a M in 
the Lorentz covariant gauge Eqn.4, we find: 

G mv (x,y) = ^ J V^Vi>i>(x)i>{y)e- $ dix ^ d ^e- w 

(16) 

Where Z is the partition function of QED4 and W is 
given by: 

e 2 f 

+ jfa)D^(x - x')j°(x') 

+ 2j^x)D^(x-x')jt(x')) (17) 

The first term in Eqn. 17 is just the conventional long- 
range four-fermion interaction mediated by the gauge 
field, it leads to the anomalous exponent in the covariant 
gauge given in Eqn. 5: 



Vi = - 



ae 
8^ 



(18) 



Note that in Landau gauge a = 0, rji simply vanishes 

The second term in Eqn. 17 is given by: 

e 2 r d A k k 2 {y - x) 2 - (1 - a)(k ■ (y - x)) 2 
(2tt) 4 J k 4 (k-(y-x)) 2 
■ (1 -cosfc- (y-x)) (19) 



Extracting the leading terms as A — > oo turns out a 
little bit more difficult than its 2 + 1 dimensional coun- 
terpart discussed in [21]. Fortunately, the a independent 
part was already given in Ref. [28] in a different context: 

e 2 3e 2 3e 2 1 

~4^ Ar + 8^ l0g Ar + 8^ (7 _ l0g2 + 2 } (20) 
Where 7 is the Eulcr constant. The linear divergence is a 
arti-fact of the momentum cut-off and should be ignored 
in Lorentz invariant regularization. 

We only need to evaluate the a dependent part: 

f A dk f 71 

ra / — / d9 sin 2 9(1 ~ cos(kr cos 9)) 
'Jo k Jo 



4^ 3 ' 



8tt 



n a l + l0g4. 

2 (logAr + 7 -2-) 



(21) 



Combining the Logarithmic terms in Eqs.20 and 21 
leads to: 

3e 2 e 2 a , . 

(22) 

Finally, we should evaluate the contribution from the 
third term in Eqn. 17 which can be written as: 



dxi/j(x)'y^ip(x) / dx' v D^ v (x - x 1 ) 



(23) 



Eqn. 23 is essentially quadratic in the fermions. Com- 
bining it with the free fermion action leads to: 

S= / dx'4>(x)(^ ll d ll + e 2 7 M / dx' v D fl „(x - x'))i[)(x) 

J J XI 

(24) 

In principal, the propagator of fermion < 
, 4>{x\)'4>{x2) > can be calculated by inverting the 
quadratic form in the above equation, but it is not easy 
to carry out in practice. Instead we can construct pertur- 
bative expansion in real space by the following Feymann 
diagrams in Fig. 3 



x' y 



(a) 



(to) 



x y 
(c) 



Fig 3: The pcrturbativc expansion series of Eq.23 

The corresponding expression is: 
G(x 1 ,x 2 ) = G {x 1 ,x 2 ) 

/rX2 
dxG (xi,x) / dx'^^D^xix - x')G a (x,x 2 ) 

/rX2 
dxGo(x!,x) dx'^^D^x- x') 

f dyG {x,y) f dy' alv D ua {y - y')G {y,x 2 ) + ■ ■ ■ (25) 
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Being just quadratic, there is no loops in the above 
Feymann series. But there may still be potential diver- 
gences. 

The explicit expression for Fig. 3b is 



F(x) 



Ai d 4 q 2 e-^ x - e 
(2tt) 4 (2tt) 4 -i( qi -q 2 )-x 
Go(qi)-l l _ l x l ,D liV (qi - q 2 )G (q2) (26) 



where x 2 — x\ = x. 

After some lengthy but straightforward manipulations, 
we can write Eqn.26 as the sum of two parts F(x) — 
F 1 {x) + F 2 {x) with 



) 4 q 2 {q~kf 



+ 



+ 



-2k ■ x7 M fc M + Pjf.Xf, 

q ■ xq 2 (q — k) 2 
k ■ - k ■ q-y^x^ 

q ■ xq 2 (q — k) 2 



(27) 



and 



4 ;„ e ikx 



1 



(2tt) 4 k 2 
d 4 q fc 2 7 M g M + g 2 7 M fc^ 
(2^) 4 (fc + <?) V 



(28) 



It is easy to see the second term in Eqn.27 is conver- 
gent. While the UV divergences in the third term in 
Eqn.27 exactly cancels as follows: 

X ■ fc7 ^- log A + X o kl /^ log A - (29) 



8tt 2 x 2 



8ir 2 x 2 



The logarithmic divergence in the first term in Eqn.27 
cancels exactly that of a independent term in Eqn.28. 
Only the a dependent divergence in Eqn.28 survives: 



ae 2 i^i^k^ 



Air 2 k 2 
The above equation leads to: 



log A 



V3 



4^2 



In all, the final anomalous exponent is 

3e 2 



V = Vl+V2+V3 



8tt 2 



(30) 



(31) 



(32) 



This is exactly the same as that calculated in the 
temporal gauge. As expected, the gauge fixing param- 
eter a drops out in the anomalous dimension r\ although 
771,772, r]3 all depend on a separately. 



C. 3 + 1 dimensional Non-Abelian SU(M) QCD 



The calculations in the last two subsections on Abelian 
QED can be straight-forwardly extended to Non-Abelian 
SU(M) QCD by paying special attention to the internal 
group structure of SU(M) group. 

In temporal gauge, the one-loop fermion self-energy 
Feymann diagram in SU(M) QCD is 




% b 



Fig 4: The fermion self-energy diagram in SU (M) QCD in tem- 
poral gauge 



The corresponding expression is 



H(K) = -ie 2 J 



Q 



(2tt)3 ll (K - Q) 2 
1 (% + ^)(T c ) ab (T c ) db 



"7j 



(33) 



where T c with c = 1, • • • , M 2 — 1 are M 2 — 1 generators of 
SU{M) group. a,d,b = 1, • • • , M are M color indices of 
fcrmions transforming as a fundamental representation 
of SU(M) group. 

It is easy to see 



E ab (QCD) = (T c y ab Z(QED) 



c\2 



(34) 



where (T c ) 2 b = C 2 (F)8 a b with the quadratic Casimir 
C 2 {F) = M 2 ^ x for the fundamental representation of 
SU(M) group. 

From Eqn.13, we have 



1 3e 2 



_ M 2 



(35) 



For SU(3) QCD, C 2 (F) = 4/3, then Vsu{3) = 
which is consistent with the result achieved previously 
with different method [3,20]. 

Eqn.35 can also be derived by using the Lorentz co- 
variant calculation presented in the last subsection. 

The results achieved in this section is not new, but 
two new, simple and effective methods are developped to 
rederive these old results. In the next section, we will 
use these new methods to derive new results in gauge 
theories with Chcrn-Simon term. 
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III. 2 + 1 DIMENSIONAL NON-COMPACT QED 
WITH CHERN-SIMON TERM 



In the two component notations suitable for describ- 
ing Time reversal and Parity breaking mass term and 
Chcrn-Simon term, the standard 2 + 1 dimensional mass- 
less Quantum Electro-Dynamics Lagrangian with Chern- 
Simon term is: 



£ = i>aln{d^ - i-j=a^)i> a + ^e^xa^ax (36) 
v N * 



Where a — 1, • • • , N are N species of Dirac fermion [26]. 

This Lagrangian was used to describe FQH transitions 
in [10,27,22,23]. As shown in [27], in straightforward per- 
turbation expansion, there is no extra UV divergences 
from the CS term in one-loop. One has to go to two 
loops to see the extra UV divergences from the CS term. 
Instead of going to two loop calculations, we resort large 
N expansion by scaling the coupling constant as g/VN 
in Eqn.36. Integrating out TV pieces of fcrmions gener- 
ates an additional dynamic quadratic term for the gauge 
field: 

+ U°(k)e^xkx)a l/ {k) (37) 
To one- loop order [27]: 

„2 



IT 



16' 



IF = 



(38) 



The dynamics of of gauge field is given by: 



+ ^{-k)U e (k)k{5^ - ^)a„(fc) (39) 

Adding the gauge fixing term ^(c^a^) 2 , we can get 
the propagator of the gauge field in the covariant gauge: 



, /,n ,£nv\kx , B knk v akuk v 



(40) 



where A = j+^.B = i^frjfe. 

Changing the last local gauge fixing term to a non-local 
gauge fixing term [25] leads to 

k tI ,k v 



/>„,/,, = - A € J^ + f (V - (1 - a)^) .41, 



In the following, we follow the same procedures in Ref. 
[21], paying special attention to the extra effects due to 
the CS term. 



A. The calculation in Temporal gauge 

With the notation K = (k ,k), in a = gauge, we 
can invert Eqn.39 to find the propagator: 



Dij(K) 



k Q 



B 
K 



«, 



(42) 



It is known that the anti-symmetric CS propagator 
does not contribute to divergence to one-loop order [27], 
the result in Ref. [21] can be directly applied by replacing 
16/Aby Sf: 



V = 



4Bg 2 
3tt 2 N 



(43) 



We expect this is the correct anomalous dimension of 
the gauge-invariant Green function Eqn.2 in 2 + 1 dimen- 
sional Chern-Simon theory. 

The tunneling density of state p(u>) ~ ui 1 ^ v . That r] 
is positive indicates the tunneling DOS increases due to 
the Chcrn-Simon interaction. 



B. Lorentz covariant calculation 

In this section, we will calculate the gauge invari- 
ant Green function directly in Lorentz covariant gauge 
Eqn.41 without resorting to the gauge dependent Green 
function. We will also compare with the result achieved 
in temporal gauge. 

First, let's see what is the contribution from the first 
term in Eqn.17. As shown in [27]. the anti-symmetric 
CS propagator in Eqn.41 does not lead to divergence, 
the symmetric part of the propagator leads to 



m = 



g 2 _B_ J_ 

N 3ir 2 



(1 - 3a/ 2) 



(44) 



In evaluating the contribution from the second term in 
Eqn.17, the extra piece due to the CS term is: 



Ag 2 



d 3 k euvxkx 



(2w) 3 N J fc 3 k 2 



'-{y - - a;) 1/(1 - cos/c • (y - x)) 

(45) 



Obviously, this extra term vanishes due to the antisym- 
metric e tensor. The result in Ref. [21] can be directly 
applied: 



m = ^(l-a/2) 



(46) 



Finally, in evaluating the contribution from Fig. 3b, the 
extra term due to the CS term is: 



F(x) c 



Ag 2 f d 3 qi d 3 q 2 e~ iqiX - e- l ^ x 
~N~ J (2tt) 3 (2tt) 3 -i( qi -q 2 )-x 
J^x v e^x(qi - <72)a 



GoGh)- 



(<7i - (iif 



G fe) (47) 
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where x^ — x\ = x. 

After straightforward manipulation, the above equa- 
tion can be simplified to: 



g 2 A r d 3 k e ikx f d 3 t 
~ l ~ J (2tt) 3 P J (2^ 



d 3 q k ■ xq 2 — k ■ qq ■ x 



) 3 q ■ xq 2 (q — k) 2 



(48) 



It is easy to find the potential UV divergences in the 
first term and second term, in fact, vanish separately ! 
Again, the result in Ref. [21] can be directly applied: 



Bg 2 a 
N -k a 

In all, the final anomalous exponent is 

ABg 2 



v=m+m+ m 



3tt 2 N 



(49) 



(50) 



This is exactly the same as that calculated in the tem- 
poral gauge. 

Although evaluating all these (1/iV) 2 corrections is be- 
yond the scope of this paper, we have established firmly 
our results Eqn.50 to order 1/N. 



C. Non-Abelian SU(M) Chern-Simon gauge theory 

It is useful to extend our results on Abelian CS 
theory established in the two previous subsections to 
Non-Abelian CS theory. Non-Abelian gauge theories 
arise both in high temperature superconductors [11] and 
FQHE [31]. In contrast to the state at filling factor at 
v = 1/2 which has a gapless Fermi surface, the state at 
filling factor 5/2 may be a gapped paired quantum Hall 
state [31]. This paired state may be a Non-Abelian state 
where the quasi-particles may obey non- Abelian statis- 
tics [31]. The effective low energy theory of non- Abelian 
states can be described by SU(2) Chern-Simon gauge 
theory. 

Extension to relativistic SU(M) Non-Abelian CS the- 
ory discussed in [30] is straightforward. By using the rule 
derived in section IIC, to the order of 1/N, we get the 
result: 



M 2 



1 



1 



7] = 



2M i + (£)2 12w 2 N 



(51) 



This results could be achieved both in temporal gauge 
and in Lorentz covariant gauge presented respectively in 
the previous two subsections. 



IV. CONCLUSIONS 

In this paper, using the simple and powerful meth- 
ods developed in [21], we calculated the gauge-invariant 
fcrmion Green function in 3 + 1 dimensional QED, 2 + 1 



dimensional QED with Chern-Simon term and their cor- 
responding SU(M) non- Abelian counter-parts by differ- 
ent methods. The calculations in temporal gauge with 
different cut-offs and Lorentz covariant calculation with 
different gauge fixing parameters a lead to the same an- 
swers. These facts strongly suggest that Eqn.13 is the 
correct exponent to one-loop and Eqn. 43 is the correct 
exponent in the leading order of 1/A^. These methods 
was previously applied to study 2 + 1 dimensional QED 
and have also been applied to many different physical sys- 
tems [29]. We summarize our results with the following 
three useful rules of thumb: 

Rule 1 : In temporal gauge, Infra-red divergence is al- 
ways in the middle of the integral and can he regularized 
by deforming the contour. The hnite exponent is the 
anomalous exponent of the gauge invariant Green func- 
tion. 

Rule 2: In Coulomb gauge, Infra-red divergence is al- 
ways at the two ends of the contour integral and can not 
be regularized by deforming the contour. The anomalous 
exponent is even not defined, but the IR divergence will 
be canceled in any gauge-invariant physical quantity. 

Rule 3: The anomalous exponent of gauge-invariant 
Green function is equal to the sum of the exponent of 
gauge-dependent Green function in Landau gauge and 
the exponent of the inserted Dirac string also in Landau 
gauge. 

Although we demonstrated the above three rules of 
thumb only to one-loop or to the order 1/A^. We expect 
they all hold to any loops or to any order of 1/N. 

I thank J. K. Jain for helpful discussions. 
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